Analysis by x-ray microtomography of a granular packing undergoing compaction 
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Several acquisitions of X-ray microtomograpliy have been performed on a beads packing while it 
compacts under vertical vibrations. An image analysis allows to study the evolution of the packing 
structure during its progressive densification. In particular, the volume distribution of the pores 
reveals a large tail, compatible to an exponential law, which slowly reduces as the system gets more 
compact. This is quite consistent, for large pores, with the free volume theory. These results are also 
in very good agreement with those obtained by a previous numerical model of granular compaction. 

PACS numbers: 45.70.Cc, 41.60.Ap 



Understanding the slow dynamics of out-of-equilibrium 
systems is still an open and debated issue. Even though 
granular media are not thermal systems, due to the in- 
significance of the thermal energy fc^T in comparison to 
the energy needed to move a macroscopic grain, they 
are of great interest in this context. Indeed, it has been 
suggested that the relaxation of a granular medium un- 
der weak mechanical perturbations, such as shaking or 
shearing, has a formal analogy with the slow dynamics of 
out-of-equilibrium thermal systems 0,13 ■ This analogy is 
based on the idea that the geometry of the system is more 
important than any other parameter as the driving en- 
ergy or the interaction between particles. Consequently, 
granular media are often presented as ideal systems for 
studying the relaxation towards equilibrium |3l|. Never- 
theless, one may keep in mind that, unlike the thermal 
energy, the mechanical agitation of the granular mate- 
rial is neither stochastic non generally isotropic and thus 
leads to different spatial scales of the energies . 

One of the most common behavior of a shaken granular 
packing is its slow compaction. This compaction is easily 
observed through the progressive increase of the volume 
fraction and appears to be a simple way to investigate 
the above analogy. The first experiments on compaction 
have been realized in a thin tube giving rise to a high 
lateral steric constraint in the packing [3 . The evolution 
of the volume fraction reveals a very slow dynamics and 
can be reasonably well fit by an empirical law in the in- 
verse of the logarithm of the number of shakes. More re- 
cently, Philippe and Bideau carried out new experiments 
with a reduced lateral confinement . They showed that 
the resulting compaction dynamics is consistent with a 
stretched exponential, also called Kholrausch- Williams- 
Watts function, which is the most usual description of 
the relaxation in a out-of-equilibrium system such as a 
glass (see for instance and references therein). 
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Contrary to most previous studies on granular com- 
paction in which the quantities measured deal with the 
properties of, at least, apart of the packing (as the mean 
packing fraction $ jj, [g or the vertical profile (§j), the 
study discussed here benefits from x-ray microtomogra- 
phy to characterize the microstructure of a packing un- 
dergoing compaction. A previous work on x-ray microto- 
mography applied to granular materials has been carried 
out recently that involves a structural study of a unique 
packing of about 2000 nearly monosized beads 8] . In this 
paper we focus on the microstructural evolution of a com- 
pacting granular bed by means of an original statistical 
analysis and we relate this to the free volume theory. 

We use an experimental setup close to previous ones 
1^ 01 : 200-400 fim diameter glass beads are poured to 
about 80 mm height in a 8 mm inner-diameter glass 
cylinder. The whole is vertically shaken by sinusoidal 
excitation at a frequency of 70 Hz. The applied acceler- 
ation is measured by an accelerometer and the intensity 
of the vibration is characterized by F, the maximal ap- 
plied acceleration normalized by gravity {g = 9.81 ms~^). 
The experiments are performed as follows: starting from 
an initial loose reproducible configuration ($ « 0.57), a 
packing is vibrated for a given number N of oscillations 
with a fixed acceleration F and then analyzed by x-ray 
microtomography. Experimental data acquisitions have 
been recorded at the ID19 beam line of the ESRF (Euro- 
pean Synchrotron Radiation Facilities). Synchrotron ra- 
diation improves upon traditional x-ray techniques and 
offers alternative methods of imaging. State-of-the-art 
microtomography provides the three-dimensional map- 
ping of the linear absorption coefficient in the bulk of 
millimeter-sized samples with a spatial resolution of the 
order of 1 /im. A description of this experimental method 
may be found in 9] . A monochromatic coherent beam is 
used to get sample radiography for 1200 angular sample 
positions ranging from 0° to 180°. An x-ray energy of 
51 keV was selected to ensure a high enough signal to 
noise ratio. The exposure time is 1 s by projection. The 
energy is selected using a classical double monocrystal 
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FIG. 1: Part of the 3D reconstruction of the initial packing 
(<E> « 0.57). 

device. A filtered backprojection algorithm 0, is 
used to compute the three-dimensional mapping of the 
linear absorption coefficient in the sample. The light de- 
tector used is based on the the fast read out low noise 
charge couple device (FRELON CCD) camera developed 
by the ESRF detector group (1024 x 1024 elements, 14 
bits dynamic). A thin scintillation layer deposited on 
glass converts x-rays to visible light. Light optics mag- 
nify the image of the scintillator and project it onto the 
CCD. With such a setup, the resolution is 9.81 by 
pixel. Each reconstructed packing contains around 15000 
grains. An example is reported in Fig^ 

Using an image processing software, the size and the 
location of each particle are recovered. Then a complete 
set of information is available to study the geometry of 
the packing. In the following we note (R) and (w) the 
mean radius and mean volume of a sphere, respectively. 

The first characteristics reported are the classical pair 
correlation function g{r). Figure|2Ia) represents the vari- 
ation of the pair correlation function before and after 
more than 10^ excitations of intensity F = 3.0. Note 
that the initial packing fraction, measured with the image 
processing software, is equal to 0.57, whereas the steady- 
state packing fraction, obtained after a large number of 
excitations, is approximately 0.62. It appears that the 
different peaks of g{r) slowly move with the number of 
taps. In particular the second peak slightly progresses 
from r — 4(_R) to r = 2\/3(i?), which is characteristic of 
developing structure. Besides, one can also remark that 
unlike Seidler et al. we do not observe a double peak 
around r = 4(7?). This is due to large size distribution 
of the grains in our experimental work (in 8] it is signifi- 
cantly weaker.) Nevertheless the evolutions observed for 
g{r) are not significant. It is even more striking if one 
compares g{r) for steady-states corresponding to differ- 
ent intensities of excitation as reported in Figl^^b). Here 
the steady-state pair correlation functions obtained for 



r = 0.95, 1.6 and 3.0 have been plotted; their packing 
fractions are, respectively, 0.617, 0.639, and 0.624. These 
correlation functions do not differ significantly. It means 
that the pair correlation function is not a suitable tool 
to study structural changes during granular compaction. 
Compaction is then a process more subtle than a simple 
decrease of the mean distance between grains. Since this 
function is averaged over all the angles of orientation, an- 
gular reorganizations cannot be seen. In contrast to the 
numerical work of Rosato and Yacoub j we do not ob- 
serve the appearance of secondary sharp peaks in the pair 
correlation function. This may be due to the difference 
in terms of vibration intensity F (up to 3 in this work, 
6.4 in 12]) and to the small size of the vessel they used, 
which can induce strong wall effects as crystallization. 
Note that by using a local bond order parameter we 
do not observe crystallization due to the walls. This is 
not surprising since the size distribution of our grains is 
rather large. 

Further evidence of a transformation in the packing 
microstructure is provided by the study of the size of 
the interstitial voids. This can also be an experimental 
test for the free volume theory 0, Q| that postulates 
an exponential decay for the distribution of these voids 
p{v) cx exp(— u/tio). Following previous works [Tsl IT^ 
we define a pore size via the Voronoi' tessellation. A 
Voronoi' polyhedron around a sphere is the region of 
space in which all points are closer to this sphere than to 
any other sphere in the packing. The Voronoi network 
that is the whole collection of the vertices and of the 
edges of the polyhedra maps the pore space. Each vertex 
of this tessellation is equidistant to four neighboring 
spheres and therefore defines the center of a pore. The 
volume of a pore is then the size of a virtual sphere 
centered on the vertex and in contact with the four 
neighboring spheres. The volume of this void-sphere 
partially refiects the volume of the whole void situated 
inside the tetrahedron formed by the centers of four 
neighboring spheres. Since in our experiments the 
spheres are not perfectly monosized, we have adapted 
this method to a polydisperse packing replacing the 
Voronoi tessellation by the so-called navigation map 
It is then possible to compute the size distribution of 
the pores for each packing. Following Philippe and 
Bideau 1^ the volume of the pores v is normalized by 
the mean volume of a grain (w) . Figure presents 
the evolution of the volume distribution of the pores 
for F — 3.0 at different stages of the compaction. 
First of all, we observe that unlike the pair correlation 
function, p{v / {w)) drastically changes with the number 
of excitations. Indeed an exponential decay law is found 
for the distribution of the voids, at least for the pores 
larger than octahedral pores (w > (^2 - If w 0.0711). 
The exponential shape persists during the compaction 
of the packing yet with a reduction of the tail, i.e. with 
a decrease of the characteristic volume Wq. A typical 
evolution of Vq is shown in the inset of Fig|3Ia). As 
reported in Fig Olb), the final exponential decays of 
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FIG. 2: (a) Evolution of the pair correlation function for F — 3.0. (b) Steady-state pair correlation functions obtained for 
F = 0.95, F = 1.6 and F = 3.0. 
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FIG. 3: (a) Evolution of the volume distribution of the pores during compaction with F = 3.0. Inset of (a): decrease of vo/{w} 
with the number of oscillations A'^ and for F = 3.0. (b) Volume distributions of the pores for the initial packing and for three 
different steady-state packings obtained for F = 0.95, 1.6 and 3.0 (b). 



the steady-state distribution, obtained after a sufficient 
number of oscillations, are strongly dependent on the 
intensity F of the vibrations. By contrast, the part of 
the distribution corresponding to the smallest pores is 
less r dependent. The remarkable point is that these 
results are very close to those obtained by numerical 
simulation by Philippe and Bideau |0. Indeed, as is 
illustrated in Fig. 01 they also found an exponential 
decay for p{v/{w)) which slightly reduces as the system 
compacts. Their simulation is solely based on the 
geometric constraints in a three-dimensional packing of 
identical hard spheres. Each tap is decomposed into a 
vertical dilation of intensity e (z z[l + e]) followed 
by a nonsequential deposition. This last stage is ruled 
by gravity (the individual motions of the particles are 
preferentially oriented along the vertical) and by the 
condition of non interpenetration between particles. 
When the gravitational energy of the whole packing 
gets nearly stabilized, the deposition is stopped and 
another tap is implemented. More details can be found 
in p^. The agreement between this simulation and our 
experiments emphasizes the fundamental importance 
of the geometry in granular compaction. The steric 



constraint between the grains seems to rule the local 
rearrangements of grains allowed by the shaking energy. 
One can observe that the small pore distributions for 
the numerical simulations and the experiments are quite 
different. We explain this point by the fact that contrary 
to our experiments the grains used for the numerical 
simulations are perfectly monosized and moreover they 
are never in real contact (see 0|.) Indeed the grain 
size distribution decreases the mean volume of the pores 
and therefore the number of small pores is larger in the 
experiments. We can also observe that the existence 
of an exponential tail does not seem to be sensitive to 
the polydispersity. However, we should point out that 
contrary to our experiments and to previous work j^, 
significant packing fraction gradients are present in 
the simulations. They are probably consequences of 
the z-dependent dilation used in the model. In order 
to quantify the importance of the tail in the pore 
size distribution we also have plotted on the inset of 
Fig. 01 vq/{w) as a function of 1/$ for experiments 
and numerical simulations, for any available number 
of taps and any shaking acceleration. We observe a 
linear correlation between the two quantities which is 
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FIG. 4; Volume distributions of the pores obtained with the 
numerical model of fl^ at different stages of a packing under- 
going compaction with the dilation parameter e = 0.1. Here 
w is the volume of the particles. Inset : plot of vo/{w) as a 
function of 1/$ for all experiments (o) and a fully represen- 
tative set of simulations (•). 



in agreement with the free volume theory 0, . This 
curve also shows that the packing fraction 4> and the 
characteristic volume vq vary in the same way during 
the compaction. 

As demonstrated here, the use of X-ray microtomog- 
raphy to study granular compaction provides much 
information on the microstructure of a granular packing 
undergoing compaction. This further opens several 
perspectives to this work. First of all, our analysis 
gives an experimental measure of the free volume and 
can be used as an experimental test of the free volume 
theory for granular compaction Then, a quan- 

titative comparison between these microtomography 
experiments and the numerical simulations described 



above can be done. In particular, we plan to correlate 
the experimental excitation parameter F with e, the 
numerical one. For that purpose we will use the measure 
of the characteristic volume of a pore vq and assume 
this quantity provides the link between F and e. We 
also plan to study the hysteresis on the steady-states 
and the microstructure of a granular material during 
the compaction along the well-known reversible and 
irreversible branches |l8| . 

In this article, we have reported experimental results 
on the development of the microstructure in a granular 
packing submitted to sinusoidal vibrations. It has been 
shown that, even at rather high acceleration (F = 3), no 
significant changes are observed in the pair correlation 
function. The structural signature of compaction is 
then more subtle than an average decrease of the mean 
distance between grains. The reorganization of the 
grains can be advantageously analyzed through the 
distribution of the pores volumes. It has been found 
on a consistent way with numerical simulations that 
compaction is mainly due to a decrease of the number 
of the largest pores. The volumes of these pores are 
statistically distributed along a broad exponential tail 
which progressively reduces while the structure gets 
more compact. 
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